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Low X physics and Structure Functions 
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Abstract. In this talk I review the behaviour of structure functions at low values of 
Bjorken x and discuss the theoretical underpinnings with particular attention to resum- 
mation schemes. I present the need for 'less inclusive' events to distinguish between various 
resummation schemes and discuss the various difficulties in differentiating experimentally 
between different schemes. 

1. Introduction 

The wealth of data on structure functions at low x that has emerged from HERA 
in the last six years or so have provided fresh and interesting challenges to theorists 
to explain the steep rise in the structure function F2{x) along with many other 
new aspects of perturbative QCD at these hitherto unprobed kinematic regions. In 
this talk I will concentrate on structure functions and discuss in some detail the 
experimental data along with the theoretical attempts made to explain them. I 
will not be touching upon many of the other interesting aspects of low x QCD viz 
Large Rapidity Gap Events, Quarkonia Production, spin physics. Photon Structure 
Functions and so on, most of which will be addressed by other speakers in this 
workshop. 

The basic interaction picture is the usual DIS diagram for the process 

e(fc) +p{p) e{k')+X. 
The following kinematic variables completely describe the process: 

s = (fc + pf ^ 4E,Ep (1) 
Q2 = _g2 _ 2E,E,,{1 + cose) (2) 

y = — ^l- —{I- cosO) 3) 
p.k 2Ee 

- (4) 



2p.q ys 

W^iq+pf^^Q^ + ys (5) 



here s is the CM energy squared, Ee{E'J) is the energy of the initial (scattered) 
electron, 9 the angle of the scattered electron in the lab frame with respect to 
the proton direction, is the negative of the momentum transfer squared to the 
proton, y is the fraction of the electron energy carried by the virtual photon in 
the proton rest frame, the Bjorken variable x is the fraction of the proton energy 
carried by the struck quark and W is the hadronic invariant mass. The masses of 
all particles are neglected. 
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Figure 1. Phase space eoverage of the F2 measurements 



At HERA, a 27 GeV electron beam collides head-on with an 820 GeV proton 
beam giving 

s ~ AEeEp - lO^GeV"^ 

which is much larger than hitherto obtained at fixed target experiments. (Since 
1994, the electron beam has been replaced by a positron beam at the same energy). 
As a result the two experiments HI and ZEUS can measure the structure function 
F2 in a completely new {x,Q^) domain (see Figure 1). 

The Born cross section for single photon exchange in DIS is given by 



27ra 



2 r 



2(1-2/) + 



l + R 



(6) 



where R is the photoproduction cross section ratio for longitudinal and transverse 
polarised photon, R = aL/cr- R has not been measured at HERA; therefore 
a QCD prescription based on parametrisation of parton densities is used. is 
directly measured from the scattered electron but x is calculated from and y 
and therefore depends on the experimental resolution of y. 



2. Measurement of F2(x, Q^) 

The data for F2 from HI and ZEUS 0, |^ is shown in figures 2, 3, and 4, along 
with some of the older fixed target data. It is clear that at fixed Q^, F2 rises 
with decreasing x down to the smallest Q^, the steepness of the rise decreasing 
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with decreasing Q^. Similarly at fixed x < 0.1 F2 rises with Q^, the rise becoming 
steeper as x decreases. In fact, as seen in figure 2, upto ~ 0.85GeT^^ the 
data favour a parametrisation based on a soft pomeron such as that suggested by 
Donnachie and Landshoff Q. For larger values of {Q^ > l.OGeV^) the usual 
QCD based parametrisations of GRV or MRS || give a good description of the 
data, as seen in Figure 3. 
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Figure 2. F2{x,Q^) as a function of x for low . The Donnachie Landshoff (DOLA) 
and the standard QCD evolution represented by GRV are shown 

The behaviour of F2 as a function of for selected values of x is presented in 
figure 4. The growth of the structure function with for x < 0.01 over three 
orders of magnitude in is clearly seen. With the new data from HERA the gap 
between fixed target and the HERA experiments has been filled and the remarkable 
agreement with a NLO QCD fit performed by HI over three decades in and over 
four decades in x is clearly seen. 

The HI experiment has parametrised the rise of F2 with a; by -F x^^ at fixed 
values. The result for A as a function of is shown in Figure 5. It is clear from 
the figure that A increases with increasing reaching a value of around 0.4 — 0.5 
around ~ 10^ - lO^Ge^^. 

3. The behaviour of F2(x) 

The general behaviour of F2{x) can be understood as follows. Below ~ 0.85GeT^^ 
the Regge theory predictions for a soft pomeron like behaviour seems to work quite 
well In this approach, the total cross section (for pp, K'^p, w^p etc. are fitted 
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Figure 3. F2{x, Q^) as a function of x for higher Q'^ compared with model predictions, 
to a form 

atot = As' + Bs-'^ (7) 

where e = 0.08 corresponds to a soft pomeron intercept a]p{Q) = 1 + 6 and r] = 0.45 
corresponds to the p meson trajectory intercept a = 1 + -q. Assuming a dipole 
model for the proton, the x dependence of F2 can be extracted; 

pDOLAr _ y. -0.08 , o 0.45 (c,\ 

Thus the DOLA analysis predicts a gentle rise of a;"" ''* for with decreasing x. 

However this prediction, as we have seen from the data fails for > 0.85GeF^ 
where standard QCD evolution takes over, given generically by the GRV fits. For 
the rest of the section we shall be discussing the nature of this 'good' agreement of 
the data with QCD evolution. 
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Figure 4. Compilation of F2 measurements as a function for selected values of x. 
The curves are a NLO QCD fit by HI. 

In QCD, the structure function ^2(2^, Q^) is given in terms of the parton distri- 
butions 



(9) 



The evolution of the structure function F2 as a function of is governed, in pertur- 
bative QCD, by the Dokshitser-Gribov-Lipatov-AltareUi-Parisi (DGLAP) equation 



ahiQ2 



2tt 



P P 

J^qq ^qg 
P P 

^99 



(10) 



Here q and g are the quark and gluon distributions, ® denotes convolution with 
respect to x, i.e.[/ ® g]{x) = ^f{-)g{y), nj is the number of flavours of quarks, 
ei is the electric charge of the quar^ qi. The coefficient functions C's at leading 
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Figure 5. Variation of the exponent A from fits of tlie form F2 ~ x' ^ at fixed values 
and X < 0.1. 



order are 



C,ix,Q^) = 5{l-x); Cg{x,Q'-)=Q. 



(11) 



At higher orders, they depend on the specific factorization scheme. In the common 
parton scheme the above remains true to all orders. 
The splitting functions have the expansion 



zvr ■' 



(12) 



The first two terms above define the NLO DGLAP evolution. 

The rise of F2 as a function of x is explained by the usual DGLAP evolution 
and also by more non-conventional dynamics like the BFKL evolution ||^. Let us 
briefly review the two different approaches. 

Perturbative QCD does not provide any information on hadronic structure. How- 
ever the idea of factorisation comes to the rescue. Very briefly, the electron proton 
cross section can be written as a convolution of the partonic cross section of the 
partons in the proton scattering off the electron times the probability of finding the 
parton in the proton. Symbolically 



(13) 
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where aei the partonic cross section is calculable in perturbative QCD. Unfortu- 
nately this cross section is plagued with infra red divergences which are of two 
kinds - soft divergences (fc 0) which are cancelled by the virtual diagrams and 
coUinear divergences (fcy — » 0) which are absorbed (factorised) into the bare parton 
distributions. The result of this last step is to get new redefined parton distri- 
butions which depend now on the factorisation scale /if. Thus the IR safe cross 
section for ep has the form 



'^•^v = 51 [/':/p(^i^) ® (14) 

■i 

The scale is arbitrary but usually taken to be equal to in which case the 
lowest order contribution to a a is the Born graph. Then fi/p{Q'^) is the parton 
density in the proton as seen by a photon of virtuality Q^. 

Since Oaiy^p) can be calculated perturbative for any scale, the change in the 
redefined parton distribution can also be calculated. Thus while perturbative QCD 
does not predict the form of /j/p, it does predict the change of fijy with scale. These 
are given by the various evolution equations, the differences between them being 
the approximations that are used to restrict the phase space for radiation. These 
approximations are then valid in regions of x, where the selected contributions 
are the dominant ones. These therefore give different evolution equations (D)GLAP, 
DLL, BFKL .... 

To understand the leading log summation technique in the DGLAP approach, 
we use the Dokshitzer method At low values of a;, the gluon is the dominant 
parton in the proton. He showed that the GLAP summation of the (a.^ logQ^)" 
terms amounts to a sum of gluon ladder diagrams of the type shown in Figure 6 (a) 
with n gluon rungs. The form of the gluon splitting function Pgg ~ Q/x for small x 
gives the evolution of the gluon distribution function as 



i-2 ■ ■ ■ / 1.2 / j,2 

^ f'l!L...fk fk^,gi^^,^Ql) (15) 

It is clear that the (logQ^)" builds up from the nested k integrations. In fact this 
contribution comes from a region where the transverse momentum of the emitted 
gluons are strongly ordered 

>> klj, » ... » klj. » kfj. (16) 

This kinematic region is therefore relevant when is large but x is not so small 
that as{Q^)ln{l/x) is large i.e 

as{Q^)Hl/x) « a,(g2)ln(QVg2) < 1 
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Figure 6. A gluon ladder diagram that contributes to the DLL or BFKL ladder sum- 
mations. 



Similarly, the (log(l/a;))" term comes from the region where the longitudinal mo- 
mentum fractions are strongly ordered 



a;«^„ «... « 6 



(17) 



In the region of small x the log(l/a;) terms have also to be summed. With the 
assumption of ^ig{^i,Qo) approaching a constant, say Go the result we get is 



Go exp < 2 



In 



Of 



TT \QiJ \X 



In 



Go } (18) 



in the limit of large and small x. 

This is an example of an all order double leading log (DLL) summation of 
as logQ^ log(l/x) obtained by summing the strongly ordered gluon ladder diagrams 
as we have just explained. It tells us that as x decreases, xg{x, Q^) increases faster 
than any power of log(l/a;) but slower than any power of 1/x. It is valid in a region 
where 



a,(Q2)ln(Q7Q2) 
a,{Q^)ln{l/x) 



« as(Q')ln(QVQo)Ml/a;) < 1 



If kiT ~ ki+iT then we lose one power of logQ^. However, at small x, this can 

be compensated by a large log(l/.T) factor. Relaxing the strong ordering constraint 
on the kx's and summing just the log(l/a;) in the small x region will give us the 
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leading log summation in log(l/a;) instead of DLL and we get 



exp 



A log 







c lo 








3)1 









(19) 



where A — (3as/7r)c. This difhcult summation was done rigorously by Balitskij, 
Fadin, Kuraev and Lipatov (BFKL) The constant c = 4 log 2 and so, for 
as « 0.2, A « 0.5. We have highly simplified the discussion of the BFKL equation, 
and in practise one needs to work in terms of the unintegrated gluon distribution 
/(x, /cy) and integrate to get the final gluon distribution. The BFKL 

region is therefore one where x is small and is not large enough to reach the 
DLL regime. 

We thus have two different predictions for the rise of F2, the DLL (from the 
DGLAP evolution equations) and the BFKL. In principle it should therefore be 
possible to distinguish between these predictions by looking at the F2 data from 
HERA. Unfortunately this is not that simple. The steepness of the rise of F2 with 
decreasing x can be controlled by varying Qq, or the starting distribution .g(^i, Ql)- 
In addition, the set of equations based only on F2 measurements is underconstrained 
and one needs one other measurement like the longitudinal structure function Fi 
in order to constrain the system fully. 

The net consequence of the above is that the present data on F2 does not dis- 
tinguish betweem the BFKL and DGLAP (or for that matter the CCFM which 
embodies both) predictions. 

One of the reasons that the data on F2 is not sufficient is due to a phenomenon 
called fcy diffusion. Let us explain this in some detail. The BFKL equation predicts 
the behaviour of the unintegrated gluon distribution f{x, k^) to be 



/(^,4) 



oc 



{x/xo 



-\/27rA' \n{x()/x) 



exp 



ln^(4/fc|) 
2 A ln(a;o/a;) 



(20) 



for fixed a,. Here A = ^41n2 ~ 0.5 for a, = 0.19, A" = ^28C(3). Thus the 
unintegrated gluon distribution exhibits a Gaussian in In k'^ with a width which 
increases with the BFKL evolution length y^\n{xo/x). This form implies that indi- 
vidual evolution paths follow a kind of randon walk in and an ensemble of paths 
exhibits a diffusion pattern according to a Gaussian in In/c^. The net consequence 
of this is that fc^ can diffuse into a region of very small values where perturbation 
theory no longer holds. Contrast this with kr ordering in DGLAP which ensures 
that the fcy's along a ladder are always contrained to remain within a region suffi- 
ciently removed from the infrared. The calculation of F2 involving an integration 
over all fcy in the BFKL formalism could therefore include diffusion of kx into the 
IR region where perturbation theory fails. A possible solution is to study special 
final state configurations where diffusion into the IR region is minimised by fixing 
the start and end point of evolution for the above region. 
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This is the idea behind a suggestion by MueUer |^ on measuring an observable 
that is less inclusive than the F2 measurement (in which none of the properties of 
the hadronic final state is measured). He suggests studying DIS events at small x 
containing an identified high-energy jet emitted close to the jet of proton fragments 
(see Figure 6(b)). The identified jet originates from the parton labelled a. If we 
study events with Xj large (> 0.05) and x very small [k, 10"'') the ratio x/xj will 
be sufficiently small to reveal the (x/xj)^^ behaviour of the BKFL ladder. Details 
may be found in the references mentioned above. This study is currently under 
way and preliminary results are reported in [p!^ . We will have more to say on this 
later. 

The behaviour of F2 as a function of both and x (for large and small x) 
can be neatly combined through the double asymptotic scaling analysis of Ball and 
Forte jl^. They build on an old analysis of de Rujula et al. jl^ to show that the 
gluon distribution function for decreasing x increases faster than a power of ln(l/a;) 
but slower than a power oil/x. There are numerous other evolution equations and 
details of the current phenomenology of these evolution equations can be found in 

One of the by-products of the measurement of the structure function F2 and its 
log evolution is that the gluon distribution inside the proton has been measured 
with greater accuracy (c/. Figure 14 in |l^ ). 

To summarise this section, the structure function F2 has been found to rise 
steeply with decreasing x down to as low as 0.85GeV^. This rise can be 
parametrised by a form x^^ where A is found to be a function of and reaches a 
value of around 0.4 — 0.5 around lOOOGeT^^. This behaviour is described very well 
by conventional NLO QCD (through the DGLAP equation). However this does 
not rule out the existence of non-conventional dynamics like BFKL because the 
system is underconstrained. Thus some other measurement which is less inclusive 
than F2 like the forward jet mentioned above, or a measurement of the longitudinal 
structure function Fi is required. In any case, some new physics must take over 
at sufficiently low x otherwise the indefinite rise of F2 would violate the unitarity 
bound. 



4. Longitudinal Structure Fiinction 

The longitudinal structure function Fl — F2 — 2xFi is related to the differential 
cross section by 



(Pa Ana^ 



dxdQ^ xQ* 



(l-y+Y)^2-^F, 



(21) 



Fl is related to F2 and the gluon distribution function through the Altarelli- 
Martinelli equation 



^F2{z, Q') + 8Y, e?(l - x/z)zg{z, Q^) 



(22) 



At low X the dominant term is the second term. 



talk February 1, 2008 9:9 



11 

A direct measurement of Fl (or R) at fixed (cc, Q^) requires data at two different 
values of y i.e. two different values of the CM energy (Q^ = sxy). Since HERA 
operates at fixed CM energy, the value of is extracted in four steps. One assumes 
that F2 is well described by NLO-QCD. Thus F2 is determined from a NLO-QCD 
fit at a low value of y (< 0.35) where the Fl contribution is negligible. This is then 
extrapolated to high values of y (~ 0.7) and then this value is used to subtract the 
F2 contribution from the cross section at ?/ ~ 0.7. A typical extraction of Fl is 
shown in the figure 7 below. 
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Figure 7. A typical extraction of Fi at y = 0.7. The shaded band corresponds to the 
allowed range in Fl from parton distributions determined from QCD as explained in the 
text. The dashed lines correspond to limits of = and -Fi = -F2. 



Before we discuss issues of more exclusive events in order to differentiate between 
various evolution equations like DGLAP and BFKL, some discussion on the details 
of the phenomenology of these equations arc in order. We have seen earlier the NLO- 
DGLAP indeed seems to work very well in describing the data down to 2GeV^. Do 
we therefore even need any other evolution equations like BFKL to describe the 
data on F2. The answer is yes, because the 'perfect' fit of DCLAP to the data 
depends crucially on the input parton densities at low Q^. To understand this 
fact, we need to compare the LO and NLO parton input densities. Thus the fact 
that the GRV parametrisation works so well is crucially dependent on the fact that 
the gluon densities have to be varied a lot. This is clear from figure 8 where the 
input gluon densities for LO and NLO are shown in the GRV parametrisation. The 
evolution of F2 is given by 

dF-, 
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Figure 8. Input gluon and u quark dcnsitcs at LO and NLO for GRV paramctrisation 



where 



■ 2.2 — ^- 

TT X 



(24) 



Notice that the LO term is fiat while the NLO term is steep in x. Thus a stable 
evolution of F-2 requires that the NLO steepness of Pqg has to be compensated by 
a gluon density which is less steep at NLO than LO. This is clear from Figure 
8 and tells us that there is a strong correlation between x-shapes of Pqg and g, 
implying that fixed order DGLAP for F2 needs extreme flexibility of parton density 
parametrisations to offset the large NLO perturbative corrections at small x. The 
NLO steepness of Pqg is in fact the lowest order manifestation of NLO BFKL in the 
quark channel. 

Recently progress has been reported in calculating the next-to-leading-log (NLL) 
contributions to the BFKL equation. Details may be found in [ p^ . Preliminary 
calculations seems to indicate that NLL contributions stem the leading order rise 
of the structure function {x^^ with A ~ 0.5) to A ~ 2.65q!,5(Q)(1 — cag{Q)) with 
c = 3.5. 
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The lesson to learn from all this is that we need detailed investigations into 
the phenomenological and conceptual issues regarding high energy log summations 
and more accurate signatures for BFKL evolution in the hadronic final state. The 
selection criteria for these have to be such that they suppress the phase space for 
DGLAP evolution and maximise it for BFKL evolution. 



5. Forward Jet Events 

One of these methods - the forward jet measurement has already been discussed. 
Preliminary results from HI are shown in Figure 9 and show a distinct preference 
for BFKL evolution, and is not reproduced by the DGLAP Monte Carlos. 
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Figure 9. Forward jet cross section at HERA as a function of x. The solid histogram 
is the NLO MEP.JET (DGLAP) result for the scale choice = At|. = ^(O.SJ^/cr)^ 
with 5=1. The two dotted histograms show the uncertainty of the NLO prediction, 
corresponding to a variation of ^ between 0.1 and 10. The BFKL result of Bartels et al 
|l7{ is shown as the dashed histogram. The data points are the new HI measurements ^] 



The situation however is still far from clear. There are few theoretical calcu- 
lations in resummed DGLAP. Most predictions for the final state observables are 
derived from Monte Carlo models whose drawback is their complexity and flexibil- 
ity to model hadronization. This makes it difficult to pin down which feature of 
the theoretical input is being tested when compared with the data for the various 
Monte Carlos in the market like MEPJET, HERWIG, LEPTO etc. On the other 
hand no Monte Carlo exist yet for BFKL evolution though some attempts have 
been made in that direction by ARIADNE. For a detailed discussion of forward jet 
events and its relevance to BFKL phenomenology, sec |16[ . 
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6. Single Particle Spectrum at large px 

Many of the ambiguities inlierent in tlie forward jet event calculations mentioned 
above can be overcome by considering single particle emission at relatively large 
Pt in the central region. Such an event is more immune to hadronization effects 
and more directly reflects In^ diffusion from the BFKL ladder. The theoretical 
calculation has been done in |1§|| and preliminary data is shown in figure 10 for one 
value of X. There seems to be some evidence of ln(l/a;) effects typical of BFKL 




Figure 10. The transverse momentum spectrum of charged particles (tt^ ,tt^ , , K^) 
in the pscudorapidity interval 0.5 < r] < 1.5 in the virtual photon-proton ccntre-of-mass 
frame. The results are shown for kinematic bin 1 with the central values x = 1.6 x 10~* 
and = 7 GcV^. The continuous and the dashed curve show the spectra obtained 
with <&,: and / calculated from the BFKL equation. They only differ in the choice of 
fragmentation scale: for the continuous curve the fragmentation functions were evaluated 
at scale ij? = (2fcj)^ and for the dashed curve at scale /i^ = When BFKL radiation 
is neglected in the calculation of the pr spectra, i.e. when the quark box approximation 
<E>,: = <I>|''' is used and strong ordering at the gluon vertex is assumed, then the dash-dotted 
and dotted curves are obtained. The fragmentation functions were evaluated at scales 2kj 
and kj respectively. The data points shown are from the HI measurement of the charged 
particle spectra 
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evolution and diffusion. However it would be useful to compare this with the full 
fixed order calculation which would allow a clearer distinction between the different 
predictions. The full DGLAP calculation for this process has not yet been done 
and is presently under study, . 



7. Conclusions 



In this talk we have concentrated on the various measurements of the structure 
funcions F2 as also Fi. We have compared the F2 measurements with various theo- 
retical calculations corresponding to different evolution equations (DGLAP, BFKL 
..). We have argued that the F2 measurements alone, being 'too inclusive' are not 
sufficient to distinguish between different resummation programs. We present var- 
ious other proposals made in the literature for 'less inclusive' events - in particular 
we discuss forward jet and single particle pr spectrum as signatures for distin- 
guishing fixed order calculations like DGLAP with BFKL and indicate that while 
there seems to be evidence of the In 1 /x type of evolution and fcy diffusion typical 
of BFKL, more theoretical calculations and more precise experimental results are 
needed before any conclusive evidence is obtained of the presence of BFKL effects. 
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